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1. Introduction 



Distribution amplitudes are wave functions of quarks and gluons within a hadron at small 
transverse separations between the constituents. They are particularly useful for the calculation of 
hard exclusive reactions at large momentum transfer. In this case, it is sufficient to consider the 
valence Fock state, i.e. the three-quark state for baryons. We have computed the normalization 
constants and moments of the distribution amplitude of both the nucleon and its negative parity 
partner, the A'^*(1535), from first principles in lattice QCD. An update on these calculations is 
presented below. 

In the continuum, in the infinite-momentum-frame, the leading-twist component of the nu- 
cleon (and similarly the A'^*(1535)) wave function can be written as 

|^,t)=/^ / -^^fiL{\uHx,y{x2)dHx,))-\uHx,)dHx2)uHx,))}. (1.1) 

J 2v24xiX2X3 

Here, the arrows indicate the spin, the transverse momentum components have been integrated out, 
Xi are the longitudinal momentum fractions, integration is defined by J[dx] = Jq dxidx2<ixj5{l — 
xi —X2 — X3), fy is the normalization constant which corresponds to the "wave function at the ori- 
gin", and is the leading-twist nucleon distribution amplitude (NDA): it contains information 
about the distribution of the momenta of the three valence quarks. 

The NDA can be expanded into (one-loop) multiplicatively renormalizable terms 

(p{xi\pL ) =120X1X2X3-^ l+Cio(xi -2X2+X3)L% +Cii(xi -X3)L''''o 

+ C20 [1 +7 (x2 -2x1X3 -2x1)] +C21 (1 -4x2) (xi -X3)L''ft 
+C22 [3 - 9x2 + 8x1 - 12x1x3] L'W + . . . I , 

where L = as{pi) / (Xs{piQ) and the Cij are the so-called "shape parameters" at a scale /Xq. In the 
asymptotic case /i — 00, all shape parameters are zero. The shape parameters of first and second 
order, c\j and C2j, are directly related to the first and second moments of the distribution amplitude, 
^/mn (-^yjtjj / + m + « = 1 or 2). The moments of the NDA are defined by 

<p'™ = |[dx]xix^x'^(p(xi,X2,X3). 

Momentum conservation (xi +X2 +X3 = 1) demands that the 9''"" fulfill the constraints 

where cp*^^^ = 1. The five shape parameters of first and second order and four constraints of the 



form eq. (1.2) uniquely determine the nine first and second order moments, and vice versa. 

The next-to-leading twist components of the nucleon wave function provide information about 
the orbital angular momentum of the quark in the nucleon. In this work, we have only determined 
the corresponding normalization constants Ai and A2. 
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Table 1: Multiplets of three-quark operators, grouped by dimension and irreducible representation 

2. Distribution Amplitudes from Lattice QCD 

To determine the NDAs from lattice QCD, several tasks have to be performed. Irreducibly 
transforming multiplets of three-quark operators have been worked out in ^ and are shown in Ta- 
ble [l[ The renormalization constants were computed in These are based on a non-perturbative 
renormalization and subsequent 1 -loop-conversion to the MS scheme. The calculation of 2-loop 
conversion factors, based on a consistent subtraction scheme for three-quark operators in dimen- 
sional regularisation, is currently in progress. 

On the lattice, correlation functions of the form 

have been calculated. Here, ^ is a smeared nucleon interpolator and ^ is a local three-quark 
operator with up to two derivatives, belongs to one of the irreducibly transforming multiplets of 
three-quark operators. The separation of positive and negative parity states, i.e. of the nucleon and 
the A'^*( 1535), has been achieved with the generalized Lee-Leinweber parity "projector" ^(1 it ^74) 

The normalization constants fy (or /n*), Ai and A2 are obtained from i^s without derivatives. 
Out of the several possible operators (see Table |I|), those which are convenient to calculate have 
been chosen. The first and second moments of the distribution amplitude require iff s with one and 

two derivatives, respectively. In order to avoid mixing with operators of lower dimension, one is 

12 4 
limited to one-derivative operators in t— and two-derivative operators in Tj- Operators with three 

or more derivatives, which would yield higher moments of the NDA, are excluded from this study. 

The first and second order shape parameters are then determined from a fit of the first and second 

moments, constrained by eq. ([T^). 

The lattices that we have used to determine the NDAs are shown in Table ^ The three volumes 

at jS = 5.29, K = 0.13632 provide good insight into finite volume effects. Also, a new lattice at a 

pion mass of less than 200 MeV has been analyzed. This brings us closer to the physical point, 

making extrapolations easier. Unfortunately, we have not analyzed lattices with reasonably low 

pion masses at lattice spacings other than a 0.072 fm (j8 = 5.29) yet. So, while the larger 

mass data of the j8 = 5.40 lattices do not indicate significant discretization effects, there is still a 

remaining continuum extrapolation uncertainty at small quark masses. 
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K 


mji/ MeV 


Size 


# Configurations 


j3 = 5.29,a = 0.0716fm 


0.13590 


660^ 


24^ X 48 


901 


0.13620 


428 


24^ X 48 


850 


0.13632 


306 


24^ X 48 


540 


0.13632 


295 


32^ X 64 


950 


0.13632 


288 


40^ X 64 


858 


0.13640 


182 


48^ X 64 


798* 


j3 =5.40, a = 0.0604 fm 


0.13610 


7221" 


24^ X 48 


687 


0.13625 


6221" 


24^ X 48 


1180 


0.13640 


503 


24^ X 48 


1037 



Table 2: Lattice parameters. The pion mass is at the given lattice volume, i.e. not extrapolated 
to infinite volume. * Results from this lattice are to be considered preliminary. ^ The pion mass 
of these lattices is beyond the range of applicability of chiral perturbation theory; therefore, these 
lattices were not used to fit the ;^PT low-energy constants. 



For the normalization constants and distribution amplitudes of the nucleon, we have results 
from 1-loop covariant baryon chiral perturbation theory to improve both the extrapolation to the 
physical pion mass and to infinite volume For each observable, there are two low energy 
constants that we fit with lattice data. For the A'^*(1535), however, results from chiral perturbation 
theory do not exist as yet. 



3. Results 

In the following, we show selected results. The leading-twist normalization constants /at and 
fy* are displayed in Figure |l^ While fy and fy* are - within the error bars - almost identical 
at higher pion masses, they are significantly different from each other at mj^ < 300 MeV: the 
nucleon's leading-twist wave function at the origin is bigger than that of its parity partner. 



The chiral perturbation theory fit to fy is shown in Figure [lb|. Once the finite volume correc- 
tions have been applied, the lattice data points lie nicely on the fitted curve. 

In Figure ^ the next-to-leading twist normalization constants, Ai and A2, of the nucleon and 
the A'^*(1535) can be seen. For the nucleon, we find that A2 ~ — 2Ai, which should hold exactly in 
the non-relativistic limit. The A'^*(1535) shows a curious behavior towards lower pion masses, for 
which we do not have a conclusive explanation as yet. 

In Figure ^ we show the volume dependence of the shape parameters. No statistically sig- 
nificant effects can be seen in most cases. One exception is C22 of the A^*(1535), where a larger 
deviation between the 32^ x 64 and 40^ x 64 lattices occurs but, given the relatively stable behavior 
of the other shape parameters, we attribute this to a larger-than-normal statistical fluctuation. 

The first-order shape parameters ci, of the N* ( 1535) are much bigger than those of the nucleon, 
indicating a stronger deviation from the asymptotic wave function. The C2i, however, (maybe with 
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Figure 1 : The leading-twist normalization constants fy and fy* as functions of the pion mass. The 
chiral extrapolation of is shown in (b). The green bands are the 1- and 2-a errors of the chiral 
extrapolation. The red points are the lattice values and the blue points are "finite volume corrected" 
values. The left most (black) data point shows the value for ff^ from QCD sum rule calculations 
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Figure 2: The next-to-leading twist normalization constants Ai (left) and A2 (right) as a function of 
the pion mass. 



the exception of the A'^*(1535) C21) are consistent with zero, within the error bars. 

Using the first order shape parameters of the j3 = 5.29, K" = 0. 13632,40^ x 64 lattice, we have 
created the barycentric plots of the wave functions in Figure ^ At asymptotically large momentum 
transfers, the wave functions will be completely symmetric. A shift of the maximum of the nucleon 
wave function towards higher values of xi (see Eq. ( |1.1| )) can be clearly seen and, in the case of the 
A^*(1535), this shift is even more pronounced. 
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Figure 3: Volume dependence of the shape parameters at j8 = 5.29 and K = 0.13632. Shown on 
the ;c-axis is the inverse length of the lattice l/L, where L = 24,32,40. 




Figure 4: Bary centric plot of the nucleon (left) and the A'^*(1535) (right) wave functions. Due to 
the large error of the second moments, only the first moments are shown in this plot. 

4. Conclusions 

The wave function normalization constants and low moments of the leading twist distribution 
amplitudes of the nucleon and A'^*(1535) have been calculated from lattice QCD. The set of lat- 
tices that were used spans a wide range of pion masses and lattice volumes. Extrapolations to the 
physical pion mass and infinite volume were made, using chiral perturbation theory formulae, in 
the case of the nucleon. 
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